Quantitative Sum Rule Analysis of Low- Temperature Spectral Functions 
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We analyze QCD and Weinberg-type sum rules in a low-temperature pion gas using vector and 
axial-vector spectral functions following from the model-independent chiral-mixing scheme. To- 
ward this end we employ recently constructed vacuum spectral functions with ground and first- 
excited states in both channels and a universal perturbative continuum; they quantitatively describe 
hadronic r-decay data and satisfy vacuum sum rules. These features facilitate the implementation of 
chiral mixing without further assumptions, and lead to in-medium spectral functions which exhibit 
a mutual tendency of compensating resonance and dip structures, suggestive for an approach toward 
structureless distributions. In the sum rule analysis, we account for pion mass corrections, which 
turn out to be significant. While the Weinberg sum rules remain satisfied even at high tempera- 
tures, the numerical evaluation of the QCD sum rules for vector and axial-vector channels reveals 
significant deviations setting in for temperatures beyond ~140 MeV, suggestive of additional physics 
beyond low-energy chiral pion dynamics. 



I. INTRODUCTION 

The Lagrangian of quantum chromodynamics (QCD) 
is invariant under the transformations described by the 
chiral symmetry group SUL{Nf) x SUn(Nf), where Nf is 
the number of light-quark flavors [D-Q ■ In vacuum and at 
low temperatures, this symmetry is spontaneously broken 
by a non-zero expectation value of the quark condensate, 
(0| qq |0). For iso-vector hadronic excitations, this mani- 
fests itself through a large splitting of the masses of the 
chiral partner mesons p(770) and ai(1260) [|[. 

Chiral symmetry is believed to be restored at high tem- 
peratures and baryon densities. As with any symmetry, 
its breaking and restoration can be described by order 
parameters. For the chiral phase transition of QCD, 
the canonical order parameter is the quark condensate. 
Lattice-QCD computations have found that the quark 
condensate "melts" with increasing temperature Q 
until its disappearance signals the restoration of chiral 
symmetry. The gradual dissolution of the quark con- 
densate is expected to be accompanied by an emerging 
degeneracy of chiral partner hadrons. Thus, the p and a\ 
masses, or more generally, the pertinent vector and axial- 
vector spectral functions, should approach each other 
with rising temperature and finally become degenerate 
when chiral symmetry is restored. 

The interplay between the temperature dependence of 
spectral functions and condensates can be analyzed using 
finite-temperature sum rules. These relate integrals over 
spectral distributions to an operator product expansion 
(OPE) which is expressed in terms of low-energy con- 
densates. The most relevant sum rules for the present 
work are the QCD sum rules (QCDSRs) @, 1 and the 
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Weinberg-type sum rules (WSRs) [QMl l| - The former are 
channel specific, i.e., they differ for the vector and axial- 
vector channels, and they contain both chirally breaking 
and symmetric operators. The latter relate moments of 
the difference between the vector and axial-vector spec- 
tral functions to chiral order parameters. Both sets of 
sum rules provide useful tools from which either the tem- 
perature dependence of the spectral functions can be 
constrained, or by which models of temperature depen- 
dence of the spectral functions and condensates can be 
tested @. 

A low-temperature medium at vanishing chemical po- 
tential can be represented by a dilute gas of thermally 
excited pions. In Rcf. [l2j], it was shown that the per- 
tinent in-medium vector and axial-vector spectral dis- 
tributions can be deduced model-independently as lin- 
ear combinations of the vacuum distributions, where the 
strength of the "mixin g" i s determined by thermal pion 
exchange (cf. also ^ was further noted that 

the mixing effect straightforwardly satisfies the WSRs, 
which has been verified subsequently [IJ [l7|. Con- 
versely, the finite-temperature QCDSRs have not re- 
ceived much attention in the context of chiral mixing. 
Rather, the (low-) temperature dependence of the con- 
densates has been evaluated and subsequently used to 
infer properties of vector and axial-vector spectral func- 
tions fla - |20| . More recently, chirally mixed spectral 
functions have been implemented within so-called finite- 
energy sum rules (FESRs) which are obtained from mo- 
ments of the QCDSRs [HUH- 

In the present work, we revisit the use of chirally mixed 
spectral functions in both WSRs and QCDSRs utilizing 
updated vacuum spectral functions (23j and numerically 
evaluating the consequences of finite pion mass correc- 
tions. Certain features of the vacuum spectral functions, 
such as their degenerate continuum in vector and ax- 
ialvector channels and quantitative agreement r-decay 
data 0, HI, as well QCDSRs and WSRs, will turn out 
to be important in a quantitative low-temeprature anal- 
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ysis. 

The remainder of this paper is organized as follows. 
In Sec. HJ we recollect vacuum and finite-temperature 
sum rules and discuss the leading temperature depen- 
dence and pion-mass corrections of the OPE coefficients. 
In Sec. IIII( we construct the finite-temperature vector 
and axial- vector spectral functions and discuss their char- 
acteristics. In Sec. IIV1 we quantitatively evaluate the 
in-medium WSRs and QCDSRs employing the chirally 
mixed spectral functions. We conclude in Sec. [V] 



II. QCD AND WEINBERG-TYPE SUM RULES 

We begin by briefly reviewing the vacuum QCDSRs 
and WSRs (Sec. Ill Aj) and their extension to finite tem- 
perature (Sees. IIIBI and III CI respectively). A more 
detailed account for the vacuum case can be found in 
the original works 04H| an d at finite temperature in 

Refs. p. ei nam. 



A. Vacuum sum rules 

QCD sum rules utilize a dispersion relation to relate 
a spectral function to an operator product expansion. A 
commonly used form (for operators with dimension less 
than 6) reads 



M 2 J s M 2 A-l 4 M 6 



(1) 



where p(s) is the spectral function of a given channel - 
Pv(s) in the vector and p~a{s) = p A {s) + f% sS(s — m„) 
in the axial- vector channel where the bar indicates the 
inclusion of the pion pole. The ellipsis designates terms 
associated with operators of higher dimension. A Borel 
transformation has traded the space- like four-momentum 
q 2 = —Q 2 for the Borel mass M 2 . The coefficients ap- 
pearing on the right-hand side involve various conden- 
sates. For each condensate, we work to leading order 
in the quark mass, m q ; that is, we consider the effects of 
each condensate but not any finite quark mass correction. 
In the vector and axial-vector channels, one has 
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The gluon and quark condensate are denoted by (^G 
and (qq), respectively, where m q is the average current 
light-quark mass; (0\) and (Of) refer to the vector and 



axial-vector four-quark condensates. Note that the only 
difference between the vector and axial- vector channels 
is due to the chirally odd C3 terms. The four-quark con- 
densates may be written in a factorized form, 



kv,a (qqY 



(3) 



where ny,A is a parameter greater than one which ac- 
counts for the deviation from vacuum-state saturation. 
In principle, the values of k need not be identical in the 
two channels as the four-quark operators are not identi- 
cal. In Ref. 23[, it was found that the QCDSRs in both 
channels can be satisfied with a single value for k. 

The vacuum WSRs are constructed from the moments 
of the difference of the vector and axial-vector spectral 
functions, denoted by Ap(s) = py(s) — Pa( s )- They read 



(WSR-0) 
(WSR-1) 
(WSR-2) 
(WSR-3) 



Ap{s) *s = lfZ(rl)-F A , (4) 
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(5) 



Ap(s) ds = -2m q (qq) = f 2 mi ,(6) 



Ap(s)sds = -2na s (Or) . (7) 



Here, Ir 2 ) is the average squared pion radius and Fa 
is the radiative pion decay constant. We have used the 
Gell-Mann-Oakes-Renner relation in WSR-2. The chi- 
ral order parameter appearing in WSR-3 is the chiral 
symmetry-breaking four-quark condensate and can be ex- 
pressed in terms of the vector and axial- vector four-quark 
condensates as 
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This condensate can also be be written in a factorized 



form as (Of B ) 



k (qq) with k 



B. Finite temperature QCD sum rules 

At finite temperatures, we will only consider mesons 
with vanishing three-momentum in the thermal rest 
frame. The finite-temperature QCDSRs then have the 
same structure as in Eq. (fTJ) but with the condensates de- 
veloping a temperature dependence and the appearance 
of additional non-scalar condensates [l8| . The latter can 
be characterized by their twist, which is defined as the 
difference of an operator's dimension and its spin. The 
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finite-temperature ci and C3 terms are then given by 



+ (O t ^ s=2 ) t , (9) 
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+ (0^= 2 ) T , (11) 



where (C T 



%, (CT= 2 < S = 4 ) T , and <cr 



are 

the twist-2 spin-2, twist-2 spin-4, and twist-4 spin-2 
condensates, respectively. A full discussion of these 
and other non-scalar condensates can be found, e.g., in 

Refs. mm. 

For a consistent implementation of the chirally mixed 
spectral functions, we need to obtain the leading-order 
temperature dependence of the condensates. To do so, 
we must introduce a small parameter about which we are 
expanding. In the current setup, there are three scales, 
m„, T, and A x , where the latter corresponds to a typi- 
cal hadronic energy scale (also representing the conden- 
sates). At low temperature and small pion mass, there 
arc thus two small ratios available, m^j A x and T/A x . 
For the temperatures in which we are interested, m n ~ T 
are of similar size. We define a parameter 



A 



m n T 

V V 



(12) 



to represent both ratios. Thereby we will work to leading 
order in A. 

The temperature dependence of a generic operator, 0, 
can be expressed in terms of its vacuum value plus its 
expectation value within external hadronic states. The 
matrix elements with singlc-pion states provide the dom- 
inant contribution at low temperatures; thus the temper- 
ature dependence can be written as 



(O)t- (0)q + 3 



d 3 k 

r(E k ;T) (n a (k)\0\n a (k)) 

(27TJ 2E k 

(13) 

where f n denotes the thermal Bose distribution and 
Eu = \J k 2 + m 2 the pion energy. Soft pion theorems can 
then be applied so that (7r a (fc)|C|7r a (fc)) becomes momen- 
tum independent for scalar operators [l8|*. Corrections 
to the soft pion theorems are suppressed by at least one 
power of A. Therefore the entire temperature dependence 



Non-scalar operators can still depend on the pion momentum 
because of their spin structure. However, the dependence is usu- 
ally simple enough that the integral over momentum can still be 
performed. 



is contained in the quantity (adopting standard notation) 

1 d 3 p 



2 f d 6 v 



(14) 



The function e(T) is proportional to the scalar pion den- 
sity, such that e oc p%j ' f 2 represents a small parameter. 
Formally it is C(A 2 ), thus the leading finite temperature 
effects to the condensates are of this order. In the chiral 
limit, = 0, it simplifies to e(T) = T 2 /Qf 2 ; however, 
for finite pion mass, we should evaluate the full expres- 
sion in Eq. (fl"4)) . 

The temperature dependence of an operator thus 
amounts to evaluating the matrix element (7r a |C|7r a ). For 
quark and gluon condensates to order A 2 it was derived 
in Ref. as 



(qq) T = {qq) [ 1 - -e(T) ) , 



(15) 



Gl) T = (^Gl) Q - \ m lfleiT). (16) 

The temperature dependence of the four-quark conden- 
sates was also found in Ref. [l8[ for k = 1. Upon express- 
ing their vacuum values with Ky,A 7^ 1, one obtains 

(0\) T = k v (qq)l (l - y^e(T)) , (17) 

18 K 
11 KA 



(o£) T = KA mi (1 



e(T) . (18) 



We note [ill, |l7J, |l8|, |27[ that factorization of the 4-quark 
condensate is not expected to hold at finite temperature; 
the appearance of Ky/A m the above expressions is a mere 
relic of the paramctrization in the vacuum. The temper- 
ature dependence of the non-scalar operators (again to 
order A 2 ) is given by E3, Ef| 



( T=2, S =2\ T = Al\-mlh{T) + I 2 {T)) , (19) 
(CT= 2 ' S = 4 > T = -A\ (jmMT) 

+ ^m 2 I 2 (T) + 2I 3 (T) 
(0^= 2 ) T = -Bl(\mlh{T)+I 2 {l 



where 



In(T) 



d 3 ( p 2)7i-l 

r(E p ;T). (22) 



2tt 3 E p 



Note that h{T) = f 2 e(T)/2. The quantities A* and A\ 
represent moments of the light-quark distribution func- 
tions, tp(x,n 2 ), inside the pion as 

An = 2 ( dx x 71 ' 1 [ip(x, n 2 ) + i>(x, n 2 )] . (23) 
Jo 
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we will use 



0.97 and A\ 



Following Refs. [lj||2? 
0.255. The parameter B^ can be determined from deep 
inelastic scattering, which in the case of the nucleon gives 
B2 = — 0.247GeV 2 . A similar measurement for the pion 
is not available. Therefore we will proceed by taking 
B2 = and assessing deviations from this value as part 
of the uncertainty in our calculations. 

Since we will need to numerically evaluate the QCD- 
SRs, a metric is needed to quantify the agreement of the 
two sides of a given QCDSR. We use the method pro- 
posed by Lcinweber [29j which constructs the average 
deviation between the left-hand-side (LHS) and right- 
hand-side (RHS) of the sum rule over a range of M 2 * 1 . 
This measure, known as the d value [29], [3(| , is defined 
by 



1 



AM 2 





LHS 








1 RHS 



<3M' 



(24) 



The Borel-mass window, within which the sum rule is 
expected to be valid, is bound by -M^ ax and M 2 lin . The 
size of the window is given by AM 2 = M^ ax — M 2 nin . 
The lower end of the window, M 2 nin , is determined such 
that the contribution of the C3 term in the QCDSR has 
reached 10% of the OPE side. The upper end, M 2 ulix , 
is found so that the continuum contribution to the LHS 
of the QCDSR has reached half of the contribution from 
the resonances. This definition is slig htly different than 
what has been used in Refs. [2!| l30j due to the larger 
continuum threshold in the spectral functions used here, 
but the resulting Borel window is comparable in size to 
those of the previous works. 



C. Finite temperature Weinberg- type Sum Rules 

The extension of the WSRs to finite temperatures was 
systematicall investigated in Ref. [ll[. Corresponding 
expressions for WSR-1 and -2 were obtained while also 
deriving both the vacuum and finite-temperature forms 
of WSR-3. We are unaware of any finite-temperature 
analogue to WSR-0. For mesons at rest in the heat- 
bath frame (i.e., for vanishing three- momentum) , the in- 
medium WSRs simplify to 



Ap(s,T) 



ds 



/ Ap(s,T)ds 
Jo 

/>oo 

/ Ap{s,T)sds 
Jo 



0, (25) 
0, (26) 
-2va s {Ol B ) T . (27) 



As above, we use the notation Ap = pv — PA to in- 
clude the pion pole in the axial- vector spectral function 1 '" . 
This is to account for the pion developing a non-trivial 
self-energy in medium. The temperature dependence of 
(Of B ) was established in Refs. as 



(Of B ) T = (Or> Q [1 - 2e(T)] 



(28) 



III. 



FINITE-TEMPERATURE SPECTRAL 
FUNCTIONS 



A main ingredient to the present work is the model- 
independent leading-order temperature result for vector 
and axial-vector spectral functions. In Ref. [l2[ it was 
shown that in a dilute pion gas a vector current interacts 
with thermal pions to produce an axial-vector current 
and vice versa. The in-medium vector and axial-vector 
spectral functions then become linear combinations of the 
vacuum spectral functions, 



Pv(s,T) 
Pa(s,T) 



e(T)} p v (s,0) 
e(T)} p A (s,0) 



e(T)p A (s,0) ,(29) 
e(T)p v (s,0) .(30) 



One can derive these relations in the same way as the 
temperature dependence of the operators were deter- 
mined through an expression analogous to Eq. (|13[) . Cor- 
rections are therefore likewise suppressed by A. The tem- 
perature dependence is also determined by the same pa- 
rameter e(T) given by Eq. (jT4")) . In the context of the 
spectral functions, e(T) takes on the role of a "mixing pa- 
rameter" between the vector and axial- vector channels. 
Extrapolating the mixing prescription to high tempera- 
tures yields degenerate spectral functions for e = 1/2. 
Since the above expressions are based on low-energy ef- 
fective theory (current algebra), it is implied that one 
should mix only the non-pcrturbativc regions of the spec- 
tral functions and leave the perturbative contributions 
unaltered. 

Since the above mixing procedure only relies on vac- 
uum spectral functions, the agreement of those with vac- 
uum sum rules is essential. Here we use the vector and 
axial- vector spectral functions of Ref. [23[, which each 
contain a ground-state resonance, an excited resonance, 
and a smooth degenerate continuum. These were con- 
structed to agree with ALEPH r-decay data [24| and fur- 
ther constrained by the WSRs. Near-perfect agreement 
was achieved for WSR-1 and -2, while WSR-0 deviated 
by -1.28% and WSR-3 by -96% (the latter is actually 
rather small on the scale of the oscillations over the non- 
perturbative integration regime). The negative sign in 
the deviations indicates the axial-vector channel's contri- 
bution to be larger. Agreement with the vacuum QCD- 
SRs was found when adjusting the values of the gluon 



t The vacuum sum rules can also be expressed in terms of pa, 
t In this definition, the continuum contribution to the spectral yielding the same structure as Eqs. I|25| l- ll27ll but with (Of B ) 



functions is moved to the RHS of each sum rule. evaluated in vacuum. 
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p v /s 

0.08 r 




(a)Vector Channel 

Pa/s 




(b)Axial- Vector Channel 



FIG. 1: Vacuum spectral functions from Ref. [23| compared 
with ALEPH r-decay data [53]. The contributions to the 
total spectral function (solid lines) are due to the ground-state 
resonance (dotted curve), excited resonance (dashed curve) 
and a universal continuum (dot-dashed curve). 
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(b)Axial- Vector Channel 

FIG. 2: (Color online) Finite-temperature vector and axial- 
vector spectral functions as following from the chiral mixing 
of the vacuum ones of Fig. [T] 



condensate, (v G ^) = °- 022 GeV4 and Kv = KA = 2 - X ' 
resulting in deviations of dy — 0.24% and d,A = 0.56% 
over the Borel windows of 0.085 GeV < M < 1.47 GeV 
for the vector and 0.089 GeV < M < 1.48 GeV for the 
axial- vector channel, respectively. 

The resulting vacuum spectral functions, displayed in 
Fig. [U contain several key features. The inclusion of 
excited states in both channels moves the continuum to 
higher energies while maintaining (or even improving) 
agreement with the r-decay data. The different masses 
of the excited vector and axial- vector states (p' and a[) 
imply chiral symmetry to be still broken in this energy re- 
gion. Consequently, the perturbativc regime commences 
at higher energies with degenerate continua in the two 
channels [23| . thus consistently delineating perturbative 
and non-pcrturbative components. 

The leading-order finite-temperature spectral func- 
tions now follow by mixing the non-perturbative parts 
of the vacuum spectral functions via Eqs. (|2"9")l and (ptU|) . 
cf. Fig. [2] Note that there is no issue of mixing con- 
tinua with different thresholds, or perturbative with non- 
perturbative parts. Chiral mixing induces a mutual flat- 
tening of the oscillatory pattern of "peaks" and "valleys" 
in the vacuum spectral functions, with the peaks in one 



channel filling in the valleys in the opposite channel: the 
valley between the p and p' is filled by the a\ peak while 
the valley between the a\ and a[ is filled by the p' peak, 
while the spectral strength in the peaks is suppressed. 



IV. IN-MEDIUM SUM-RULE ANALYSIS 

In this section, we implement the finite-temperature 
spectral functions derived from chiral mixing into the 
WSRs (Sec HVAT) and QCDS Rs (Scc. UVBl) . followed by 
a brief discussion (Sec. IIV C|) . 



A. Weinberg Sum Rules 

The in- medium WSRs, Eqs. $25§-$ZT$, depend on the 
difference of the vector and axial-vector spectral func- 
tions, Ap. Using chiral mixing, Eqs. (|29l) and ([50)1 . one 
has 

Ap(s,T) = p v (s,T) - p A (s,T) 

= [l-2e(T)]Ap( fl ,0). 1 ' 
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Thus the in- medium WSR-1 and -2 remain intact pro- 
vided their vacuum limit is satisfied (in fact, residual 
deviations in the vacuum are suppressed toward chiral 
restoration). Since the leading-temperature dependence 
of the four-quark condensate is given by (1 — 2e), WSR-3 
is also satisfied by the chirally mixed spectral functions. 
This result holds for any vacuum spectral function pro- 
vided that no further adjustments are applied, e.g., in 
the treatment of the continua [Ojl . 



B. QCD Sum Rules 

Let us begin by revisiting the QCDSRs in the massless- 
pion limit of the chiral-mixing scenario. One has e = 
T 2 /6/ 2 so that the leading temperature dependence is 
of order T 2 . Both sides of the sum rule can then be 
expanded to this order. The LHS takes the form 



1 f Pv(s,0) s/M 2 

M 2 J s 

, , T s 1 f Pa(s,0) -py(s,0) _ s/M 2 
[ >M 2 



(32) 



The second integral can be evaluated using the vacuum 
QCDSRs. Because it reflects the difference between the 
axial- vector and vector spectral functions, the only term 
that survives is the chirally breaking 4-quark condensate; 
thus 



LHS 



1 

M 2 



Pv(s,0) 



-s/M 2 



^a s (OD • (33) 



For the RHS (OPE side), the vanishing pion mass makes 
the quark-condensate term and the temperature correc- 
tion to the gluon condensate vanish, while the non-scalar 
terms are higher powers of T 2 compared with e(T). The 
terms that remain are the Co term, the gluon condensate, 
and the vector 4-quark condensate along with its tem- 
perature dependence. Upon expanding the C3 term, the 
temperature corrections to the vector 4-quark condensate 
can be written in terms of the chirally breaking 4-quark 
condensate, 



■56 



-7T(X 



<oX) T 
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81 ' * /T 81 
Thus the RHS takes the form 



Tra s (Ol)+e(T)ira s (Oi B ). 

(34) 



RHS = c + 
1 



1 



' c 2 
24M* \ tt 

cr \ (35) 



Comparing the LHS and RHS of the QCDSRs one now 
sees that the chirally breaking 4-quark condensate can- 
cels between the two sides while the remaining terms are 
simply the vector QCDSR in vacuum. A similar deriva- 
tion can be made for the axial-vector channel as well. 



Thus, for masslcss pions, the QCDSRs are analytically 
satisfied to order e (see also Ref. 2l|). 

However, for finite pion mass, some of the simplifica- 
tions on the OPE side no longer occur. The expressions 
for the LHS remain unchanged, while the temperature 
dependence of the vector 4-quark condensate still can be 
expressed in turns of the chirally breaking 4-quark con- 
densate so that the cancelation between the two sides 
for this term occurs. However, terms which vanished or 
were temperature suppressed are now present, involving 
factors of the squared pion mass. Most notably, there is 
an incomplete cancelation between the temperature cor- 
rections of the quark condensate, the gluon condensate, 
and the twist-2, spin-2 condensates. A similar lack of 
cancelation applies to the dimension-6 non-scalar opera- 
tors. The balance of these partial cancelations produces 
terms on the RHS which are order A 2 where there are 
no such terms on the LHS. Thus the QCDSR are ex- 
plicitly violated at order A 2 . Therefore we must proceed 
by numerically evaluating the QCDSRs to determine in 
how far quantitative deviations affect the validity of the 
mixing scheme at a given temperature. 

The two sides of each channel's sum rule are displayed 
in Fig. [3] for four different temperatures as a function of 
Borel mass, while Tab. U quantifies the average devia- 
tion, Eq. (|24[) . over the pertinent in-medium Borel win- 
dow. Since neither channel's vacuum spectral functions 
exactly satisfy the vacuum QCDSRs, this will add resid- 
ual deviations for the in medium analysis. Therefore it 
is important to compare the in medium deviations with 
the vacuum's deviation. Overall, the d values stay within 
a factor two of the vacuum values for temperatures up 
to 120-140 MeV for the vector and axial- vector channel, 
respectively. These temperatures also roughly reflect the 
range within which the d values are approximately linear 
in e, as one would expect in the regime of validity. Thus 
one may conclude that the QCDSRs are reasonably well 
satisfied up to these temperatures. 

Let us address some of the uncertainties in the eval- 
uation of the QCDSRs. Closer inspection of the Borel 
window reveals its lower end to expand significantly with 
temperature, see Fig. [3] The decrease of A^^ in (T) can 
be attributed to the decrease of the C3 term with T, i.e., 
the 10%-level contribution of the 1/Ai 6 moves to lower 
M 2 . In fact, the low-end portion of the Borel window 
(extending outside the vacuum window) causes most of 
the increase of the deviation at higher temperature. If 
instead one postulates the in-medium Borel window to 
be the same as in vacuum (indicated by the vertical lines 
in Fig. [3]), the d values decrease appreciably at higher 
temperatures (quoted in parenthesis in Tab.|T|, thus stay- 
ing within a 1% margin up to temperatures as high as 
160 MeV. 

Further uncertainties are due to the non-scalar terms 
in the in-medium OPE. The accuracy of the twist-2 con- 
tributions, A% and A\, is not well known, while has 
not been measured for the pion. We find the impact of 
the twist-2 terms on the d values to be small, while -BJ 
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T (MeV) 





100 


110 


120 


130 


140 


150 


160 


170 


180 


e 

d v (%) 
d A (%) 



0.24 
0.56 


0.06 
0.32(0.29) 
0.65(0.57) 


0.08 

0.38(0.33) 
0.70(0.58) 


0.10 
0.48(0.39) 
0.78(0.61) 


0.13 
0.64(0.51) 
0.90(0.67) 


0.16 
0.85(0.64) 
1.08(0.76) 


0.20 
1.11(0.74) 
1.30(0.88) 


0.23 
1.43(0.97) 
1.60(1.01) 


0.28 
1.82(1.17) 
1.98(1.17) 


0.32 
2.29(1.39) 
2.53(1.34) 



TABLE I: Average deviation of the QCDSRs over the Borel window for the vector and axial- vector channels at select temper- 
atures. Values in parentheses are based on a frozen Borel window identical to the vacuum one. 



Vacuum T=100MeV T=140MeV T=180MeV 

QCDSR QCDSR QCD SR QCD SR 




(a)Vector Channel 




(b)Axial- Vector Channel 



FIG. 3: Comparison of the LHS (solid curve) and RHS (dashed curve) of the QCDSRs the vector (upper panels) and axial- vector 
(lower panels) channels at select temperatures. The extent of each plot corresponds to the Borel window at that temperature, 
while the vertical lines designate the Borel window in vacuum. 



produces a more significant deviation. For example, a 
20% change in the value of causes a relative change 
of less than 0.7% in dy and <1a for all temperatures. A 
similar variation in A\ causes a larger relative deviation 
in dy and dA of ~ 5% at the higher temperatures, but 
less than 0.1% at the lowest temperature. A larger uncer- 
tainty arises from the numerical value of -BJ. While this 
parameter has not been experimentally measured for the 
pion, a relative estimate of its size may be made by adopt- 
ing the value for the nucleon, = — 0.247GeV 2 . We 
then find that dy is altered by less than 1% at 100 MeV 
but as much as 25% at 180 MeV. The effect on the axial- 
vector channel is not as pronounced with the maximum 
relative deviation of 5% at the highest temperature. 



implementation as given by chiral mixing, due to over- 
laps of the non-perturbative resonance regime with the 
perturbative continuum in the two channels. This be- 
comes particularly problematic in the WSRs, whose low- 
temperature versions require a strict mixing of the entire 
vacuum spectral functions. In the present work, we re- 
solve this issue by employing vacuum spectral functions 
with a universal separation of non-perturbative and per- 
turbative regimes, thanks to excited states and a degen- 
erate continuum. This ensures a straightforward fulfill- 
ment of the low-temperature WSRs (i.e., within chiral 
mixing) but also maintains a numerical agreement of the 
QCDSRs in both vector and axial- vector channels at low 
temperature. 



C. Discussion 

Let us finally put our results into context with previous 
works on in-medium sum rules. Most QCDSR analyses 
thus far have been based on a ground-state plus contin- 
uum ansatz for the vector and/or axial- vector spectral 
function. This usually implied that (a) the continua in 
the two channels are different, and (b) a reduction of 
the in-medium threshold was deduced from the QCD- 
SRs. Both properties can potentially cause complica- 
tions when carried out within a rigorous low-tempcraturc 



The WSRs and the QCDSRs for chirally mixed axial- 
/vector spectral functions are analytically satisfied in the 
chiral limit, without an apparent limitation due to tem- 
perature for their validity. However, for finite pion mass, 
the QCDSR are violated at order A 2 . The then required 
numerical evaluation of the QCDSRs does exhibit devia- 
tions suggesting their breakdown at temperatures of the 
order of the pion mass. It is interesting that this break- 
down roughly sets in at temperatures where one expects 
corrections from higher resonances to become important. 
This is not seen if only massless pions are considered. 
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V. CONCLUSIONS 

In the present work we have performed a simultane- 
ous analysis of QCD and Weinberg-type sum rules with 
vector and axial-vector spectral functions in the low- 
temperature limit. Specifically, we have utilized updated 
vacuum spectral functions as an ingredient to their chi- 
ral mixing, and numerically evaluated the role of pion 
mass corrections in the sum rules. The in-medium spec- 
tral functions confirm earlier findings of a mutual flatten- 
ing trend with increasing temperature. While the WSRs 
could be shown to remain valid analytically (thereby 
avoiding a mixing of perturbative and non-perturbative 
components in the spectral functions), the finite pion 
mass corrections required a numerical analysis of the 
QCDSRs. They were found to be reasonably satisfied 
up to temperatures of about 140 MeV, with deviations 
markedly increasing beyond. This may be an indirect 
indication for additional physics, beyond low-energy chi- 
ral pion dynamics, which is, of course, well known from 



other contexts (resonance excitations). Our work may 
be extended to include more realistic medium effects, ei- 
ther by working to higher orders in T / K x and m n / A x , 
or by considering spectral functions which are known to 
agree with dilepton emission data. It is instructive to 
note that even the most simple, yet rigorous, construc- 
tion of low-temperature vector and axial- vector spectral 
functions satisfies sum rules rooted in QCD while sug- 
gesting a mechanism that could play a role toward their 
degeneration. 
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